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1. Introduction 

It is a well known fact that in the large N limit the vacnnm expectation 
value of a hnite product of observables is the product of vacuum expectation 
values. It is very easy, almost trivial, to prove this remarkable factorization 
in the Feynman diagram approach|l|. It is also possible to prove it in the 
canonical formalism [^. But in the path integral formalism it is not clear how 
this factorization holds. This factorization suggests that only one conhgura- 
tion enter into the path integral. This is true for vector models but in matrix 
models the number of degrees of freedom is order and the expansion pa¬ 
rameter of the large N limit is 1/iV^, therefore the quantum corrections to 
any classical conhguration must be order one and fluctuations to all orders 
must enter into the path integral. This explains why is so difficult to hnd 
the master held: because it is not given by one conhguration but arises from 
several non equivalent conhgurations. For instance, in one hermitian matrix 
models the partition function is given by 

Z = J D^exp{-Ntr[V{^)]}. (1) 

The potential is a function of N degrees of freedom only: the eigenvalues of 
$. Actually 

Z = J DUDAA'^exp{-Ntr[V{A)]} (2) 

where A is the diagonal matrix of the eigenvalues of d>, A is a Van der 
Monde determinant and the integration over the set of unitary matrices DU 
is trivial. The leading contribution to the partition function is given by only 
one conhguration of eigenvalues, the master held, which is a solution of the 
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saddle-point condition 

^ (ln{A2exp{-iVtr[l^(A)]}}) = 0 (3) 

In this case the factorization of observables holds trivially. 

But if we add to the matrix potential a coupling with an external matrix 
A then the potential becomes a function of degrees of freedom. Actually 

Z = jDUDAA^exp{-mrlV(A) + U^AUA]} (4) 

We can prove the factorization of a product of observables which are functions 
of the eigenvalues, like < Tr<hTr<h >. After the exact integration over DU 
the partition function becomes 

Z(A) oc J DAAexp{-Ntr[V{A) + AA]} (5) 

Only one conhguration of eigenvalues yields the leading contribution to the 
partition function in the large N limit. Therefore we can also dehne a master 
held in this case but hrst we must perform the integration over N'^ degrees 
of freedom exactly, in other words all conhgurations of the unitary matrix 
variable U enter into the master held. 

The existence of a master held for the eigenvalues explains the factoriza¬ 
tion of observables like < Tr<hTr<h >. But for observables like 

< Tr($A)Tr($A) >=< Tr{U^AUA)Tr{U^AUA) > (6) 

the factorization is difficult to understand. Of course we know that the fac¬ 
torization holds, we can prove it perturbatively, but we do not understand the 
factorization of observables nonperturbatively in the path integral approach. 
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We can expect that in the large N limit the partition function is given by 
the integral representation of hnite N but restricted over some subset of 
conhgurations. For instance (^) in the large N limit should be 

Z=[ DX [ DUA‘^exp{-Ntr[V{A)+ U^AUA]} (7) 

Jrjj 

where Fa and F[/ are subsets of the set of all eigenvalues and of all unitary 
matrices respectively. We know that there is only one conhguration in the 
subset Fa- The factorization 

< Tr^ATr^A >=< Tr^A >< Tr^A > (8) 

suggests that there is only one conhguration in Fc/. But at hrst sight this is 
not possible because the set of unitary matrices is labeled by degrees of 
freedom. But if there are several nonequivalent conhgurations in F^/ the fac¬ 
torization of observables like (^) becomes very unnatural. This old problem 
has been pointed out by Itzykson and Zuber in|^ several years ago. 

In order to understand this problem we must hnd the conhgurations which 
enter into the path integral in the large N limit. In this paper we will study 
one hermitian matrix models coupled to an external matrix. Actually we will 
study: 

Z = J D^exp{-NTr[V{^) + A^]} 

Z = J D^exp{-NTr[V{^) + A^A^]} (9) 

The hrst model has been solved in^, |[ and the second in|0. But we can 
not read from the solutions of the above models the conhgurations that yield 
the leading contribution to the large N limit. In the usual approaches one 
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must perform an integration over a subset of variables exactly. In order to 
understand the factorization of observables we must study the large N limit 
before performing any integration. 

We will show that, in the above matrix models, the configurations which 
enter into the path integral in the large N limit can be labeled with only N 
variables. For instance the subset F^/ is a subset of dimension N. Therefore 
the quantum corrections to any classical conhguration are subleading in the 
large N limit. Hence only one configuration yields the leading contribution 
to the path integral in the large N limit and the factorization of observables 
follows trivially. 

2. Integration over the unitary matrix in the large N limit 

Let us start with the following matrix model 

Z{A) = J D^exp{-N[Tr{^A)+TrV{^)]} . (10) 

Without lost of generality we can consider diagonal external matrices. 

Let us perform the change of variables 

$ ^ U^AU, ( 11 ) 

where A is the diagonal matrix of the eigenvalues of <h and 17 is a unitary 
matrix. The partition function becomes: 

Z{A) = J DAA‘^exp{-NTrV{A)}Z{A,A), (12) 

where 

Z{A,A) = J DUexp{-NTr{U^AUA)}. (13) 

This integral can be solved with the Itzykson-Zuber formula[^. But we will 
show that in the large N limit the above integral with A hxed can be reduced 
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to an integral over nnitary matrices with only order N nonvanishing elements. 
Hence, as in vector models, we can use saddle-point methods to calculate it. 

In the introduction I have claimed that at hrst sight an integral like (^) 
cannot be solved with a saddle-point method. This is because any correction 
must be order the number of degrees of freedom which, in matrix models, is 
the same as the inverse of the expansion parameter. Therefore any correction 
must be order one. But the above argument does not take into account the 
compactness of the unitary set. Unitary matrices must satisfy the constraint 

UU^ = I (14) 

Therefore, the absolute values of the elements of an arbitrary unitary matrix 
must decrease when N increase. Hence the range of variation of U decrease 
with N. But there are exceptions to this rule. For instance the absolute 
values of the non zero elements of a diagonal matrix are order one no matter 
how bigger is N. Actually for bidiagonal, tridiagonal and in general matrices 
with a hnite number of nonzero elements in each row and column, the abso¬ 
lute values of its non zero elements are order one. Therefore we can expect 
that these conhgurations must play some special role in the large N limit. 
Actually I will show that these conhgurations yield the leading contribution 
to the partition function in the large N limit. 

There are two kind of unitary matrices. Unitary matrices with order 

— N elements so smaller that they decoupled from the matrix potential in 
the large N limit, let us call them unitary matrices with N degrees of freedom, 
and unitary matrices where order N‘^ elements will be coupled between them 
through the matrix potential in the large N limit, let us call them unitary 
matrices with N‘^ degrees of freedom. Hence the partition function should 
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be the sum of two terms 


^ ^uncoupled ^coupled ( 1 ^) 

where Zuncoupied is the partition function restricted over the subset of unitary- 
matrices with N degrees of freedom and Zcoupied is restricted over the subset 
of unitary matrices with degrees of freedom. 

It is not possible to solve Zcoupied bnt it is trivial that 

Zcoupled ^ Volex^{-NTr{UlAUsA)} ( 16 ) 

where Vol is the volume of the subset of unitary matrices with degrees 
of freedom and Ug is the absolnte minimnm of the potential in the snbset of 
unitary matrices with degrees of freedom. 

The partition fnnction Zuncoupied is given by 

Zuncoupied = Aexp{-N‘^F} ( 17 ) 

where A is the factor which arises after the integration over the — N 
degrees of freedom which deconpled from the potential and I will show that 
F is the free energy of a matrix model restricted over nnitary matrices with 
only N nonzero elements. I will show that A and Vol are of the same order 
in the large N limit: this is because the compactness of the nnitary gronp. 
Then if 

F < j^Tr{U^AUsA) ( 18 ) 

then Zuncoupied gives the leading contribntion to Z. And because Zuncoupied 
can be rednced to an integral over N degrees of freedom only; actnally I will 
show that 

Zuncoupied J DUexp{-NTr{U^AUA)} ( 19 ) 
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where DU is the measure over unitary matrices with only N nonzero el¬ 
ements; there is a master held in the variables U and corrections to this 
master held will be sub leading in the large N limit. 

In order to classify the unitary matrices according to the absolute values 
of its elements we must transform the measure over the holomorphic elements 

DU = l[dUab (20) 

ab 

into a measure which includes the absolute values of the elements 


DU=\{d\Uab\ 

a<b 


( 21 ) 


The invariant measure DU can be extended over the set of all matrices 
as follows: 


JdU = J n l^dU,bdU:,6(j2UacU;,-Sa)j 

This measure is invariant under the unitary transformation. 
Let us perform the following change of variables: 


( 22 ) 


(23) 


where 6ab are the phases of the elements of the unitary matrix. The Jacobian 
of is 

N 

(24) 


N 

n 

a,6=1 


Therefore the measure (^) becomes 

N ( / N 


J DU = J f[ I d\Uab\ dOab E \Ubc\ " tdbc) " 

a,6=1 I \c=l / 


( 25 ) 
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and the partition fnnction ([I3|) becomes: 


z(aa) 


, N N 

n id\Uab\ dOab) n \Uab\ 

a,6=1 a,6=1 


N 


a,6=1 


N 


^ |h^ac||h^6c| G^p(^^ac ^^6c) 

.c=l 

N 

exp <j -N \Uab\‘^^bbAa 


dab 


a,6=1 


(26) 


The variables \Uab\ are restricted by the delta fnnctions which yield the 
following constraints: 


N 

l^acl \Ubc\ exp {iOac “ i9bc) = dab- (27) 

C=1 

Let ns stndy now the behavionr of the measnre (^5]) in the large N limit 
for differents snbset of nnitary matrices. For homogeneons conhgnrations: all 
\Uab\ are of the same order, the nnitary constraints (pT]) restrict the variables 
{|t7a6|} to be order 1 /a//V. 

Let Zh the partition fnnction (p^ ) restricted over the set of homogeneons 
conhgnrations. Then there is a trivial bonnd to Zh given by 


Zh < Volexp [-NTriUjAUsA)} , (28) 

where Ug is the absolnte minimnm of Tr{UlAUsA) in the snbset of homo¬ 
geneous conhguration and Vol is the volume of the subset of homogeneous 
conhgnrations. Because the dimension of the subset of homogeneous conhg- 
urations is the volume must be order 

Vol oc (^)^' 

WN 
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( 29 ) 






Let us study now the partition function Zi restricted over the following 
inhomogeneous conhgurations: for each row and column there are a hnite 
number of elements with absolute value order one, let us call them large 
elements, and the others are order less than 1/a/]V, let us call them small. 
Order less than Ij^fN means that the small elements are order 

|C/*r"l OC (30) 

where e is a positive number which goes to zero in the large N limit. 

The are three contribution to the constraints (pTl) . The contribution 
which comes from the large elements 

i: itinidr'i exp (ie^ - iOt,). (31) 

c 

This is order one because the sum over c can only have a hnite number of 
terms. Let us remark that these inhomogeneous matrices must have only a 
hnite number of large elements in each row a column, otherwise the above 
contribution becomes order greater than one for a = b and the matrices are 
not unitary. 

The contribution which comes from the small elements 

i: vr“"l vr““l exp - Wt,). (32) 

c 

This is order 



because now the sums over c runs over order N elements. And the contribu¬ 
tion which comes from the crossing of large and small elements. 

i: \u‘:rn\utr°“\ exp (>9„o - ^ 0 ^,). (34) 
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This is order 


1 


1 




[Vn] 


(35) 


because the sum over c is over a hnite number of terms. There are only a 
hnite number of nonzero large elements for fixed a. 

Now let us dehne e as 


e = 


Vh^ 

Then e goes to zero in the large N limit and 

f 1 




0 


(36) 


(37) 


Hence, in the large N limit the unitary constraints becomes 


E \u‘"n\ulTn exp (i9„ - ie^) = S^. (38) 

c 

Therefore the small elements have decoupled from the unitary constraints. 
Let us remark that the large elements can be put on a matrix, say Ul, and 
the constraints (^) are equivalent to the unitarity oi Ul- 

But we cannot perform the large N limit inside the Dirac delta function, 
in other words is not true that in the large N limit 


S(UU* - /) ^ 6(UlUI - I) 


(39) 


This is because the matrix Ul, the matrix of the large elements only, has only 
order N variables and there are N'^ delta functions in the partition function 

Let us consider the unitary constraints up to the hrst correction in the 
large N limit. Then 

E \U‘r’\\Ulr‘\ exp(*«„ - i»»o) +E |C'r“"l|C7r“"l<ixp - Wfc) = 5.^, 

c c 

(40) 
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The cross term (Q) gives the second correction. Because there are only N 
variables a hnite number in each row and column, the large variables 

only appear in order N constraints. For instance a matrix with N large 
variables in only one row gives contribution to all the constraints but such 
matrices are not unitary. 

Therefore, there are N delta functions which depend on the large variables 
and in the large N limit they become. 



If we call In the set given by the hrst N integers then T is a subset of the set 
In X In and has only order N elements. In other words (^) is the product 
of order N delta functions. 

And there are N'^ — N delta functions which are independent of the large 
variables: 


n dEit'a 

(a,b)Gr \ c 


'||C/£““"|exp(!9„-!9fc) 


(42) 


where f is the complementary of T in In x In- Let us remark that there 
is a complete decoupling between the large and the small variables in the 
measure. 

Let us dehne 6U as the matrix given by the small elements only. Then 
the delta functions (p!]) and (|i2|) means that 


{UL)ac{Ul)cb = (a, 6) G T 

{5U)ac{5U^)cb = 0 (a,6)ef (43) 


Let us remark that the second equation in (|43|) does not mean that 5U 
is singular because the indices a and b do not run over all integers between 
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one and N. Actually we can write the second equation as 


SU5U^ = J 


( 44 ) 


where J is an arbitrary matrix with nonzero elements Jah for the labels (a, b) 
that belongs to F and 

Jab oc (a, h) e r. (45) 


Let us consider the matrix potential now. The matrix potential can be 
split into two parts: 

N N 

y = E \Uab?KbAaa-NY. ^ (46) 

a=l bGlarge a=l bGsmall 

The hrst term is order because the sum is extended over the variables 
which are order one and the sum over the index b has only a hnite number 
of terms. In the second term each sum over indices a and b is order N, but 
the order of the absolute values are less than 1/^/N and therefore the second 
term is order less than Therefore the matrix potential is a function of 
the large variables only in the large N limit. 

Let us perform the integration over the small variables. We must take 
into account that 


5{5U5U^-J) = 


det U 
1 


5{6U^ - J{5U)-^) 
--5{5W - J{5U)-^) 


- -, (47) 

Vdet J 

Hence, taken into account the decoupling of the small variables, equations 


and (|30|) , and the dehnition of e, the integration over the small variables 
becomes 

, , Cl" Ar2_ 


1 


{{N'^-N){l+e)-eN) 
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Hence after the integration over the small variables, Zi becomes in the 
large iV limit 

Z, oc ( ^ j I DUlDUIS{UIUl - 1) exp {-NTriUlAULA)} (49) 

and after the integration over DUj^ 

Z, oc ( I DU^exp {-NTriUlAULA)} (50) 

where now DUl is the measure over unitary matrices with only order N el¬ 
ements nonzero. Let us remark that that the number of nonzero elements 
is not hxed but it is order N. Therefore the dimension of the set of inho¬ 
mogeneous matrices is not a hxed number. This means that DUl is not a 
Riemanian measure in the usual sense. This is not important in the large N 
limit. Actually I will show that the integral over DUl is given by only one 
conhguration in the large N limit. 

But we can also dehne DUl explicitly as a sum 

DUl = DiUl + D 2 UL + ■ ■ ■ (51) 

where DiUl is the measure over unitary diagonal matrices, D 2 UL is the 
measure over unitary matrices with two diagonals a so on. For instance let 
us consider the set of inhomogeneous matrices with large N elements in the 
main diagonal only. Then we perform all the tricks describe above and we 
arrive to the partition function restricted over inhomogeneous matrices with 
large elements in the main diagonal 

Z] oc ( I D.ULexp {-NTriUlAULA)} (52) 
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where DiUl ^s the measure over diagonal unitary matrices. In this way we can 
write the partition function Zi as a sum over partition function restricted over 
diagonal, bidiagonal and so on. The factors in front of each partition function 
Z] are differents, but all are of the same order in the large N limit. These 
factors are not important because in the large N limit only one conhgurations 
enter into the path integral Zi. 

The partition function Zj can be solved with a saddle-point method be¬ 
cause depends only on order N degrees of freedom. Actually the saddle-point 
equation is given by 

-^TriU^KUA) = 0 (53) 

oUab 

In this case any correction to the saddle-point solutions are order the ex¬ 
pansion parameter, which is 1/A^^, times the number of degrees of freedom, 
which are order N. Let us remark that if we perform any change of variables 
in (^), for instance 

{UL)ab^mL)ablOab). (54) 

the saddle-point conditions (|^) do not change. This is because the Jacobian 
of any change of variables grows as an exponential of N, because there are 
only N degrees of freedom, for instance the Jacobian of (|5^) is the product 
of N elements 

n \(UlU, (55) 

(a,fe)er 

whereas the potential grows as the exponential of 

Taken into account that UlU^ = /, the saddle point conhguration of (^) 
are given by: 

UAU^A = AUAU\ (56) 
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or in other words the matrix $ must commute with the external matrix A. 
Actually ( [5B|) is the condition of commutativity of <h and A in the basis where 
the matrix A is diagonal. 

Let us remark that the above equation (^) comes from the fact that Ul 
is a unitary matrix with only order N nonzero elements. It is not important 
if all its nonzero elements are order one or I/a/IV. Of course a unitary matrix 
with only order N nonzero elements must have N large elements. 

The same analysis can be carried out for conhgurations with large ele¬ 
ments, small elements and order l/y/N elements. For instance conhgurations 
with order N large elements and order — N elements order 1/'/N. Let 
us remark that there is a very important difference between elements order 
l/\fN and small elements, which are order less than l/\fN. The former are 
coupled with the large elements and between them through the potential and 
the small elements do not enter into the potential and do not coupled with 
the large elements through the unitary constraint. Nevertheless we must be 
careful with the above division between large, small and homogeneous ele¬ 
ments. An element order l/y/N enter into the potential only if there are 
order homogeneous elements in the same row or column. 

Therefore we can construct conhgurations with order degrees of free¬ 
dom coupled between them starting with the homogeneous conhgurations 
and changing a hnite number of homogeneous elements in each row an col¬ 
umn by small elements. In this way we can also construct conhgurations with 
order N degrees of freedom coupled between them starting from the inhomo¬ 
geneous conhgurations and changing a hnite number of small elements in each 
row and column by homogeneous conhgurations. Therefore with this trick 
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one can generate the same number of coupled and uncoupled configurations. 

We can also generate new configurations with degrees of freedom 
starting from the homogeneous conhgurations and adding order N large el¬ 
ements. For instance a unitary matrix with N large elements and — N 
homogeneous elements. But with this trick we can also generate inhomoge¬ 
neous conhgurations starting from inhomogeneous conhgurations: if we add 
order N large elements to a matrix with order N large elements the result 
will be a matrix with order N large elements. 

We can generate new conhgurations with order degrees of freedom 
from homogeneous conhgurations changing a complete row or column of ho¬ 
mogeneous elements by small elements. If we apply this trick to an inhomo¬ 
geneous conhguration then the new conhguration will have N large elements 
and N homogeneous elements and these will be coupled to the large elements. 
But the number of degrees of freedom will be also order N. 

All the conhgurations generated from homogeneous and inhomogeneous 
conhgurations with combinations of the above tricks exhaust the set of uni¬ 
tary matrices. Therefore for each conhguration with degrees of freedom 
generate from some homogeneous conhgurations we can generate a conhg¬ 
uration with N degrees of freedom from an inhomogeneous conhguration. 
Hence the partition function in the large N limit can be split into two terms: 

Z uncoupled T ^coupled 

where Zuncoupied is the partition function restricted over the set of unitary 
matrices with only order N nonzero elements. Actually 

Zuncoupied Ndegexp{-NTr{U^AUA} (58) 
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where ?7 is a solution of ( ^61) and N^eg is the number of solutions of (|^). 
The partition function Zcoupied is bounded by 


Zcoupled ^ Vol exp{-NTr{UlAUsA} 


(59) 


where Us is the absolute minimum in the subset of unitary matrices with N'^ 
degrees of freedom. And Vol is of the same order as Ai in the large N limit. 

Hence if the absolute minimum of the potential which is given by (^) 
belongs to the subset of conhguration with N degrees of freedom then the 
partition function Zuncoupied gives the leading contribution to the partition 
function in the large N limit. 


Let us study now the solutions of (|56|) . There are several cases: If the 
spectrum of the matrix A is non degenerate, then the matrix f/^Af/ must be a 
diagonal matrix because only a diagonal matrix can commute with a diagonal 
nondegenerate matrix. And if the spectrum of A is non degenerate, then the 
identity matrix is the only solutions of (|5UD. In this case the value of the 
matrix potential in ( [5D| ) is Tr{AA). If the spectrum of A is degenerate, then 
there are an inhnite number of unitary matrices which satisfy U'^AU = A: 
all the unitary matrices corresponding to change of basis that leave invariant 
the subspace of eigenvectors of <h with the same eigenvalue. These unitary 
matrices are block diagonal, the dimension of each block is given by the 
dimension of each subspace. Then if these dimensions are order N then 
there are solutions of (^61) that do not belong to the subset of inhomogeneous 
conhgurations and the reduction of the number of degrees of freedom does not 
take place. But this does not hold because the Van der Monde determinant, 
which arise in the change of variables (^3]), forbids conhgurations where the 
matrix A is degenerate. In the large N limit the degeneration of A cannot 
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be order N and the solutions of belong to the subset of inhomogeneous 
conhgurations. 

If the spectrum of A is degenerate, then the matrix U^AU must be block 
diagonal. For instance if 

An Ui t 1, * * * , TTi 

An = 02 i = ni, • • • , 77,2 


An 7 Tij _ ]^, * * *, TV 


( 60 ) 


then WAU is given by j boxes in the diagonal of dimensions given by the 
set of numbers Uk- In this case the matrix U is also block diagonal and 
the dimension of each block depends on both matrices A and A. If the 
degeneration of A is order N then there are solutions of (^) which do not 
belong to the subset of inhomogeneous conhgurations, but we can chose the 
external matrix A to be nondegenerate. 

In the large N limit the eigenvalue conhgurations are given by continuum 
functions a{x) and A(a:). Therefore near a given eigenvalue there are an inh- 
nite number of eigenvalues, but this number is order less than N. Therefore 
the solutions of (^) are unitary matrices which belong to the subset of inho¬ 
mogeneous conhgurations in the large N limit and the leading contribution 
to the partition function is given by Zj. 

All the solutions of (^) give the same contribution to the partition func¬ 


tion. Actually the value of the matrix potential for solutions of (p6D is always 
Tr{AA). And one expect that the number of solutions must be a continuum 
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functional of the eigenvalue configurations. Hence, in the large N limit the 
partition function is 


A^ 5 [A(a:), a(a;)] exp{Tr(AH)} 


( 61 ) 


where Ns is the number of solutions of (^). It is not possible to hnd the 
functional Ns from the equation (^61). But the integration in (|^ can be per¬ 
formed exactly for hnite iVQ. Therefore, from the exact solution is possible 
to extract the functional Ns in the large N limit. 


Let us remark a few things about the equation (56). In the application of 
one hermitian matrix models to condensed matter, mesoscopic and nuclear 
physics the matrix is identify with a Hamiltonian and the partition function 
is a sum over Hamiltonians which satisfy some symmetry [^]. For instance 
matrix models with real matrices are related to models with time-reversal and 
spin-rotation symmetries, complex matrices are related to models with time- 
reversal symmetry broken by a magnetic held or magnetic impurities, and 
quaternionic matrices are related to models with the spin-rotation symmetry 
broken by strong spin-orbit scattering for instance. In all this cases the 
correlation between the eigenvalues {Aj} are given by 


n(Ai - A,)" ( 62 ) 

where (3 can be 1, 2 or 4 depending on the symmetry of the model. Actually 
(3 is the number of degrees of freedom in the matrix elements: 1 for real 
matrices, 2 for complex matrices and 4 for quaternionic matrices. But it is 
well known that if one introduce into the matrix potential a coupling with an 
external matrix, as in ([l0|) , the exponent in (|62D changes. But let us observe 
that, in the large N limit, the coupling with an external matrix restricts 
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the sum over all hermitian matrices to the sum over hermitian matrices that 
commute with the external matrix A (|5B|). Therefore the change in the index 
(3 is still related with a change in the symmetry of the set of matrices over 
which the path integral is dehned. 

Let us perform the calculation of 

Z = J D^exp{-NTr{^A^A + V{^))} (63) 

Now the integral over the unitary matrices U is given by: 

Z = J DUexp[-NTr{U^AUAU^AUA)} (64) 

One can prove that the classical conhgurations of the matrix potential 
are given by 

{U^AUAf = {AU^AUf. (65) 

The solutions of (^) are also solutions of the above equations. But there 
are other solutions. For instance, the conhgurations which satisfy 

U^AUA = -AU^AU (66) 


are also solutions of the new 
of solutions depends on the 
is the identity matrix there 
solutions of (|BBD only if: 


saddle-point equation. In this case the number 
diagonal matrix A and A. For instance, if A 


are not solutions of (|66|). Actually, there are 


^jj 


(67) 
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for some i and j. Therefore the form of the matrix A must be: 


/ ^1 






( 68 ) 


V 

where Ai are matrices proportional to the identity matrix and D is a matrix 
which does not verify (0- Then the form of the matrix f/lAf/ must be 


/ 0 Si \ 

b\ 0 

0 B 2 

Bl 0 


(69) 


V 0 

where Bi is a complex matrix, and the square root of the eigenvalues of the 
hermitian matrices B\Bi are the eigenvalues of A. The value of the matrix 
potential for these conhgurations is 


- 2 ^ alTr{B\B,) ^ Tr{A^K^), (70) 

i 

where Ai = ail. Whereas for the same matrix A and A, the value of the 
action for the solution of (^) is 

2j2a,Tr{BlBi) = Tr{A^A^). (71) 
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Therefore, solutions of 


U^KUA = AU^KU 

U^AUA = -AU^AU (72) 


are saddle-point configurations of the partition function But their con¬ 
tributions to the partition function are different. Therefore, we can think the 
two set of solutions as different vacua of the model. There are other saddle- 
point conhgurations: conhgurations for which the matrix WAUA has two 
blocks in the diagonal, in such a way that one of the block satisfy the com¬ 
mutation conditions and the other the anticommutation conditions. These 
sets of solutions give different contribution to the integration over the U 
matrix. 

This is an interesting difference between the linear coupling to an external 
matrix and the gaussian coupling: the former has only one vacuum whereas 
in the gaussian coupling there are several vacua of the saddle point equa¬ 
tion corresponding to the integration of the unitary matrix. But in all this 
cases if the degeneration of the external matrix is order one in the large N 
limit then the solutions of saddle point equation are given by inhomogeneous 
conhgurations and the reduction of degrees of freedom holds. If the external 
matrix A does not verify (^) then the only solutions are given by solutions 
of 

U^AUA = AU^AU (73) 


In this case the eigenvalues must be all of the same sign because if for some i 
and j, the corresponding eigenvalues verify Aj = —Xj then there are solutions 
different from (??. This is because the symmetry between A and A in the 
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equations. For this vacuum the degeneration factor iV 5 [A, A] must the same 
as in the linear case. Therefore the partition function (^), in the large N 
limit, becomes: 

Z«^exp{-AfTr(A=A=)} (74) 

Even though this model can be used to dehne pure Quantum Gravity 
with extrinsic curvature its solution is not very interesting from a physical 
point of view|^, Q. 

4. Conclusions 

The main conclusion of this paper is that for matrix models with 
degrees of freedom for hnite N there is a reduction of the number of degrees 
of freedom in the large N limit, actually the number of degrees of freedom 
becomes order N. This reduction is the origin of the factorization of the 
vacuum expectation values of products of observables in the path integral 
approach in the matrix models studied in this paper. 

The reduction showed in this paper is different from the well known 
Eguchi-Kaway reduction|T^. The EK reduction is a reduction of the ex¬ 
ternal space: the model reduces to a zero dimensional model; while in the 
reduction studied in this paper the reduction takes place in the internal space 
of degrees of freedom: actually the matrix models studied in this paper are 
zero dimensional. 

It seems that this reduction is more fundamental than the EK reduction. 
This reduction is the origin of the factorization of product of observables and 
the factorization is the origin of the EK reduction. 

This result holds only if at least order degrees of freedom can be 
thought as the elements of some unitary matrices and the absolute minimum 
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of the matrix potential, as a function of this unitary matrices, is given by a 
conhguration with only order N degrees of freedom. The hrst condition is 
accomplished by every matrix model and also by gauge theories in the lattice 
because every hermitian or unitary matrix can be split into the diagonal 
matrix of its eigenvalues and a unitary matrix:<h = The second 

condition seems to depend on the matrix potential. For instance, in two 
unitary matrix models with matrix potential given by 

TrV{U,) + TrV{U 2 ) + Tr{UiU 2 ) + h.c., (75) 

there are saddle point conhgurations which do not verify the second condi¬ 
tion. But if one perform the following change of variables 

Ui ^ (76) 

where di is the diagonal matrix of the eigenvalues of Ui and is a unitary 
matrix; then the saddle point conhgurations corresponding to the matrix 
variable 17 verify the second condition. Actually this problem is analogous 
to the matrix models studied in this paper. Therefore if for any given matrix 
model or gauge theory in the lattice is possible to hnd a change of variables, 
as for instance (l76|), such that order degrees of freedom are the elements 
of some unitary matrices, for instance 17 in ([761), and perhaps there are other 
N degrees of freedom, for instance the diagonal matrix 9 in (j76|) , and the 
saddle point conhgurations of the potentials as functions of the former unitary 
matrices are given by N degrees of freedom only, then the reduction of the 
number of degrees of freedom holds in the large N limit. 

This reduction must be a general result of the large N limit because the 
factorization of the expectation values of products of observables is a general 
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result and is difficult to understand how the product of observables depending 
on degrees of freedom can factorize if fluctuation to all orders enter into 
the path integral. Therefore we can conjecture that this reduction in the 
internal space must be a general behaviour of the large N limit in the path 
integral approach. 
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